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Abstract 

The link between the nucleon generalized parton distributions and 
the non-diagonal one-body density matrix in momentum space is stud- 
ied. Attention is focussed on the region where quark generalized parton 
distributions (GPD's) describe emission and reabsorption of a single ac- 
tive quark by the target nucleon. The correct covariant connection with 
wave functions used in any constituent quark model is established. Re- 
sults obtained with different constituent quark models are presented for 
the unpolarized quark GPD's. 
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1 Introduction 

In recent years much theoretical activity has been devoted to generalized par- 
ton distributions (GPD's). They are defined as nondiagonal hadronic matrix 
elements of bilocal products of the light-front quark and gluon field operators 
and they interpolate between the inclusive physics of parton distributions and 
the exclusive limit of electromagnetic form factors |[ ||, ^, ||, ^. As such 
they contain information on parton correlations and the internal spin structure 
of the nucleon. But the main interest to study GPD's comes from the fact that, 
due to the factorization theorem 0, ||] , they are candidate to provide us with 
a unifying theoretical background suitable to describe a variety of inclusive and 
exclusive processes in the deep inelastic scattering (DIS) regime. In particular, 
they enter the cross section for exclusive photon production, i.e. deeply virtual 
Compton scattering (DVCS) ||, ||, and hard diffractive electroproduction 
of (longitudinal) vector and pseudo-scalar meson jl^. In addition, in the for- 
ward limit they become diagonal matrix elements giving the usual DIS parton 
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distribution, and their first moment gives the nucleon elastic form factors Q. 
Basic properties of the GPD's have been reviewed in Refs. |3[ 

Much effort is underway related to the measurement of these functions. A 
general discussion about the problems arising in the experimental study of the 
GPD's and the first experimental evidences have been presented in Ref. p5{ . 
Prospects for future measurements of DVCS are studied in Ref. Q. On the 
other hand, model calculations are a necessary prerequisite to define suitable 
strategies to extract GPD's from experiments. In the literature there are two 
approaches used to model the nucleon GPD's. One is a phenomenological con- 
struction based on reduction formulas where GPD's are related to the usual par- 
ton distributions by factorizing the momentum transfer dependence due to the 
nucleon electroweak form factors |l^, [l^ . This leads to double distribution 
functions parametrizing the nonforward matrix elements involved in DVCS and 
hard exclusive electroproduction processes discussed, e.g., in Ref. Another 
approach is based on direct calculation of GPD's in specific dynamical models. 
The first model calculations were performed using the MIT bag model pO[|. Fur- 
ther calculations have been performed in the chiral quark-soliton modelpl], ^ . 

A complete and exact overlap representation of GPD's has been recently 
worked out within the framework of light-cone quantisation [ p3| , . In particu- 
lar, quark GPD's are obtained as overlaps of light-cone wave functions (LCWF's) 
indicating that they are deeply connected with the non-diagonal one-body den- 
sity matrix in momentum space occurring in any nonrelativistic many-body 
problem. A preliminar investigation of such an approach was presented in 
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In fact, direct calculation of LCWF's from first principles is a difficult task. 
On the other hand, constituent quark models (CQM's) have been quite suc- 
cessful in describing the spectrum of hadrons and their low-energy dynamics. 
Therefore it is interesting to explore the connection between GPD's and CQM 
wave functions at least in the allowed kinematic range where only quark degrees 
of freedom are effective. This corresponds to the kinematic region where GPD's 
describe how a quark is taken out from the proton and, having undergone a 
hard scattering, is inserted back as a quark inside the scattered proton. This 
connection has been studied in Ref. where a simple nonrelativistic quark 
model has been used to calculate one of the twist- two GPD's at a low-energy 
scale. However, as in the case of parton distributions calculated from CQM wave 
functions a nonrelativistic approach is obviously insufficient |^ and a conse- 
quence of it is that the obtained GPD's are not defined on their natural support 
and the particle number and momentum sum rules are not always preserved. 
The support violation occurs also with the MIT bag model |20 because the 
initial and final nucleons are not good momentum eigenstates. These are com- 
mon problems arising also in the case of parton distributions calculated from 
CQM wave functions |^, ^ where they were cured by appropriately taking 
into account translation invariance | |30| . 

CQM's rely on quantum theory with a finite number of degrees of freedom. 
In this case relativity can be incorporated quite naturally by utilizing relativistic 
Hamiltonian dynamics l31| and the Bakamjian- Thomas p2| construction of the 
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Hamiltonian for a system of interacting particles. In such an approach CQM 
wave functions can be considered as eigenfunctions of the nucleon Hamiltonian 
in the instant-form dynamics and can simply be related to wave functions in 
any form of relativistic Hamiltonian dynamics. Examples are the light-front 
description of electromagnetic form factors (see Ref. 33 and references therein) 
and the covariant calculation of the nucleon electroweak form factors in 
a chiral CQM |35| with the point-form approach. Similarly, one can obtain a 
link between CQM wave functions and LCWF's with the corresponding trans- 
formation from the instant-form to the front-form representation as it is here 
proposed. GPD's in the allowed kinematic region are then obtained in a covari- 
ant approach and they exhibit the exact forward limit reproducing the parton 
distribution with the correct support and automatically fulfilling the particle 
number and momentum sum rules. 

The paper is organized as follows. In Sect. 2 the relevant definitions and 
properties of GPD's in terms of LCWF's are recalled. In Sect. 3 the connection 
is established between the (light-cone) front-form and the (canonical) instant- 
form description. This connection is used in Sect. 4 to obtain the unpolarized 
GPD's and their forward limits, i.e. the usual DIS parton distributions in terms 
of CQM wave functions. The results obtained with different CQM wave func- 
tions are presented in Sect. 5, confining some technical details in the Appendix, 
and concluding remarks are given in Sect. 6. 



2 The unpolarized generalized parton distribu- 
tions 

For definiteness let us consider virtual Compton scattering where a lepton ex- 
changes a virtual photon of momentum q^^ with a nucleon of momentum , 
producing a real photon of momentum q'^ and a recoil nucleon of momen- 
tum P'^ . Ultimately we are here interested in kinematic conditions similar to 
those familiar in the DIS regime that is characterized by the Bjorken limit, i.e. 

= —q^ — > oo, P • (7 ^ oo, and fixed xb = /2{P ■ q). When focusing 
on the deeply virtual kinematic regime of q^, i.e. DVCS, a generalization of 
the Bjorken kinematics is considered, namely the cm. energy s = [P -[- g)^ 
and the photon virtuality are large while the invariant momentum square 
t = = 2P ■ A is small ^, Under these conditions the factorization 
theorem |^, || tells us that the amplitude factorizes in a hard scattering part 
(which is exactly calculable in perturbative QCD) and a soft, nonperturbative 
nucleon structure part. The contribution of the hard scattering part to leading 
order corresponds to the so called handbag diagram with photon scattering on 
a single parton (quark or antiquark). Consequently, the soft part is a quark- 
quark correlation function, representing the process where a parton is taken out 
of the initial nucleon and reinserted back into the final nucleon after hard scat- 
tering. To describe such a process, according to Ref. Q it is useful to choose 
a symmetric frame of reference where the virtual photon momentum q^ and 
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the average nucleon momentum P'^ = ^{P^ + P"^) are collinear along the z 
axis and in opposite directions (Fig. |l]). It is also useful to use the component 
notation — [a^ ,a±] for any four- vector with light-cone components 

= {aP ± a'^)/\/2 and the transverse part di__ = (a^, a^). 





P' 



Figure 1: The symmetric frame of reference. 

In the following we only consider the case of unpolarized quarks inside the 
nucleon. Then for each flavor q the soft amplitude in Fig. |^ becomes 



, (1) 



where n is a lightlike vector proportional to (1,0,0,-1), A (A') is the hclicity 
of the initial (final) nucleon and the quark-quark correlation function is inte- 
grated along the light-cone distance at equal light-cone time (y+ — 0) and 
at zero transverse separation (y_L = 0) between the quarks. The resulting one- 
dimensional Fourier integral along the light-cone distance y~ is with respect to 
the quark light-cone momentum k = xP . The so called skewedness parameter 
^ describes the longitudinal change of the nucleon momentum, 2^ — — A+ZP"*". 
The link operator normally needed to make the definition (||) gauge invariant 
does not appear because we also choose the gauge = which reduces the 
link operator to unity. 

Following Ref. jj] the leading twist (twist-two) part of this amplitude can 
be parametrized as 

Fl^{x,i,t) = ^u(p\X')j+u{P,X)H^{x,i,t) 

+ 4+ uiP\ A') ^(P^ A) E'^ix, ^, t), (2) 

2P+ 2M 



where u(P, A) is the nucleon Dirac spinor and H'^{x, t) and E'^ix, ^, t) are the 
chiral-even (helicity conserving) and chiral-odd (helicity flip) GPD's for partons 
of flavor q, respectively. 

An explicit expression in term of LCWF's can be obtained following the lines 



4 



of Refs. The hadronic state is written as 

Mfc±]jV«'A,Ar,/3(r)|iV,/3; fci, . . . , /CA,), (3) 

N 

where the integration measures are defined by 

['^fc^^W = (n Ifefa) 2(2-)^ 6 (^f: fc,. - P^^ . (5) 

The function 'i'\^N,i3{r) is the momentum LCWF of the A^-parton Fock state 
|A'', I3;ki,. . . , /cat) (also including gluons), with /? indicating flavour, helicity and 
colour quantum numbers. The LCWF's depend on the momentum coordinates 
(collectively indicated by r) of the partons relative to the hadron momentum. 
The 7V-parton states are normalized as 

{N',p';kl...,k'j,\N,p;h,...,kN) 

N 

= 6N'NSp,f3Y[2k+{2nfS{k'^+ - k+)5{k'^, - k^i). (6) 

i=l 



j\7" a ' 



f 


dx 







Correspondingly, the hadron states are covariantly normalized as 

(P', A'|P, A) = 2P+(27r)3 5{P'+ - P+) 5{P'^ - P^) 6x'x, (7) 

with 

J2 fmN[dk±]N\'^x,NAr)\^ = 1- (8) 

Having in mind the link between GPD's and CQM wave functions we have 
to confine our discussion to the region ^ < a; < 1. In this region and in the 
symmetric frame 

N,I3 j=l 



X / [cfx]Ar[dfcj_]jv5(5- Xj)*A',A',/3('"')*A,7V,/3(r)e(Xj),(9) 



where Sj labels the quantum numbers of the jth parton, P specifies all other 
quantum necessary for the iV-parton state, and the set of kinematic variables 
r, r' are defined as follows: for the final struck quark, 
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for the final TV — 1 spectators {i ^ j), 



1-e 

and for the initial struck quark 



1-C 



1 — X.i 



P 



for the initial iV — 1 spectators (i 7^ j), 



Working out the spinor products we have 



:£;«(s,e,i), 



(11) 



(12) 



(13) 



(14) 
(15) 



where 



and 



Ai + iA^ 
|Ax| ' 

4^2 M2 



(16) 



(17) 



- 1-^2 

is the minimal value for — i at given ^. 

Therefore, one can extract H'^ and E"^ separately from the knowledge of 
Fy^. In particular, E'^ is directly given by Eq. (p^), and 



iJ'?(x,e,t) 



F^+(x,?,t) + 



(18) 



3 Nucleon wave functions in front and instant 
form 

In this section a connection will be established between the (light-cone) front- 
form and the (canonical) instant-form description. 

When necessary, labels [/] and [c] on the wave function will refer to front 
form and canonical form, respectively, and vectors with tilde are defined in the 
front form. 
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3.1 Single-parton states 

Omitting colour degrees of freedom that do not matter in the following consid- 
erations, single-parton states can be defined either with front-form coordinates, 

IM,t)[^], fc = (fc+,fci), (19) 
(with isospin r) or with instant-form coordinates, 

|fc,A,r)[^], k = {k^,ky,k^), (20) 
with corresponding normalizations 

[/](fc',A',r'|fc,A,r)[^] = 2k+ {27Tf6{k+ ~ k'+) S{h - k'^) Sxx' 6^^' , 

[,](fc',A',T'|fc,A,T)[^] = 6{k-k')SxySrr'- (21) 

The transformation from the instant- to the front-form representation reads 
IM,t)[^] =^ / rffc'|fc',A',T')[,][c](fc',A',r'|fc,A,r)[/], (22) 

A'r' 

where, according to Eq. (4.35) of Ref. Q, 

[,](fc',A',T'|fc,A,T)[/] =\/2^(27r)3/2 5(fc-fc')i^i('(i?c/(fc))5rr'. (23) 

Here Rcf is a Melosh rotation, and the normalization constant \/2a;(27r)'^/^ with 
uj = ko = (fc"*" + k~)/V2 derives from a different normalization used in Ref. | ]3l[ . 
Then Eq. ( p2| ) becomes 

|fc,A,r)[^] = V2^(27r)3/2^i?^/^(i?,/(fc))|fc,A',T)[,]. (24) 

A' 

3.2 N-parton states in the front form 

In the front-form description the iV-parton contribution \N;P^X)y^ to the nu- 
cleon wave function can be derived from the zero-momentum state by a suitable 
boost Lf{Q) with Q = \P\/Mo {Mq being the mass of the noninteracting N- 
parton system) and 

Lf{QUl,0,0,Or = —. (25) 
In fact, under a Lorentz transformation A we have 

C/(A)|7V;P,A)[^] -^|Af;AP,A')[y](A'|i?v^(A,P)|A), (26) 

A' 
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where i?vi^(A, P) is the Wigner rotation associated with the A transformation. 
Assuming for A the boost A = LJ^{Q), the Wigner rotation becomes the iden- 
tity, Rw{Lj^iQ),P) = 1- Thus 



U{Lfm\N;0,X)[f]^\N;P,X) 



'[/]' 



U{Lj\Qj)\N;P,X)^^^^\N;0,X)^fy 
In turn, according to Eq. (pi). 



|iV;0,A)[^] = V2Mo(2^)'/'|iV;0,A)[,]. 
On the other hand, according to Eq. (||) we have 



(27) 
(28) 

(29) 



/ 


dx 







N 



N 



dki_]N'^^^\{x^, fc±,; X,,n}) Yl |fc„ A„ r,)[^.j. (30) 



3.3 Connection with the instant-form representation 

In order to estabhsh the connection with the instant-form representation one 
has to work out the deha functions and the measures (Eqs. ^ and (^)). Since 
J2i = Mq, we have 



N 



N 



N 



In addition, 

N 



N 



i=l 



akzi Ui 



(31) 
(32) 

(33) 



We can then rewrite the zero-momentum state d3(J) as 



|iV;0,A) 



[/] 



E 



N 



Y[dk, 



X,, 1 1 



LU, 2(27r)3 



N 



2{2nfMoS 



N 



X«'[fl({x„ A„T,}) II V2^(2^)3/2 
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E 



TV 



1/2' 



N 



Mod [J2h] ^^^\{x,,k^f,K,n}) 



N 



;2(27r)3 n[2(2^)3]-i/2 ^ ^1/2^ ^„ r,) f^,. (34) 



4=1 



Therefore, the zero-momentum state in the canonical representation, 

i^;0'^)[c] = E / Il'^^'M E^O *A'({fc-A-^a)ni^-^-^*)[c]' (35) 

ri,Ai i=l \i=l / 1=1 

can be rewritten in the hght-front representation using the definition ( p9| ) and 
the resuh fpl) as 



|iV;0, A), 



E 



N 



- \ 1/2' 



^5 (^E^') *[fl({x„fci,;A„Tj) 



AT 



: [2(27r)=^] V2 [][2(2^)3]-V2 ^ Z?i/2^(i?,^(fc,))|^,, r,)^^,. (36) 



From Eqs. ( |3q ) and ( pq ) a relationship can be obtained between hadron wave 
functions ^'l^' and This relation is here worked out in the case of a nucleon 
with only three valence quarks, i.e. iV = 3. The three-quark (valence) Fock state 
of the nucleon has been shown in Ref. to have only one independent LCWF 
for all configurations where the quark helicities add up to the helicity of the 
nucleon. This state is also known to almost exhaust the contribution to GPD's 
at large values of x ||2^ . 

By definition we have 



5' (fci , fc2 , fcs ; Ml , n , M2 , T2 , M3 , T3) 

= (fci,/ii,Ti; fc2,/i2,T2; fc3,M3,T3|^;0, A)[c]. 



(37) 



The bracket on the r.h.s. can be expressed in terms of the light-front represen- 
tation (p6|), so that 



*A (^l'^2,fc3;Ail,Tl,/^2,T2,Ai3,T3) 



1 



2(27r) 



X1X2X3 



UJ1UJ2UJ3 



1/2 



X E ^ifM(i?c/(^l))i?;./L(i?c/(^2))i?if,3(i?c/(fc3)) 
Ai A2 A3 

x'i'\'{xi,k±i,X2, k±2,X3, fcj.3; Ai, Ti, A2, r2, A3, T3). 



(38) 
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Viceversa, 



\ f] -* 

W^'{xi,k±i,X2,k±2,X3,k±3;Xl,Ti,X2,T2,X3,T3) 



2(27r) 



1 aJiW2'^3 



Mo a:iX2a;3 



1/2 



xD;,g;(i?,,(fc3))*l:i(^i,fc2,fc3 ;Ml,n,Ai2,T2,M3,T3)- (39) 



Eqs. ( p8[ ) and (|39|) give the correct covariant transformation linking wave 
functions obtained in the (canonical) instant form, e.g. in some CQM, and the 
corresponding valence-quark component of the LCWF. 



4 The valence-quark contribution 

The valence-quark contribution to GPD's is obtained by speciahzing Eq. (^ to 
the case iV = 3, i.e. 



^A'A(3^,e,i) 



xvI/W*(/,{AJ,{r,})vI/yi(r,{AJ,{n})e(x,), (40) 



where the front-form wave function ^Pl; (r, {A;}, {r^}) is related to the corre- 
sponding wave function \l/[^'({fci}, {A.;}, {ri}) in the canonical form by Eq. (^. 
Separating the spin-isospin from the space part of the canonical wave function, 



'f^i{k^}, {A,}, {n}) = ^(fci, fc2, fc3)$Ar(Al, A2, A3, n,T2, T3) , (41) 



we have 



f[^l(r,{A,},M) = 2(2^)3 



1 LUlUJ2l-^3 



1/2 



ip{ki,k2,k3) 



Mo X1X2X3 

E <i:(^c/(fci))<i*(^^/(fc2))<i:(i?c/(fc3)) 



X $Ar(A*l,M2,Ai3,Tl,T2,T3), 



(42) 



where the Melosh rotations are given by 

Dli^Rcfik)) - {X\R,f{xMo,kM 



, . m -t- xMq — i<7 ■ {zx k±) 

iM ; — IM/- 

{m + xMo)^ + kl 



(43) 
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4.1 Parton distributions 

In the limit — )■ 0, where x goes over to the parton momentum fraction x, some 
of the GPD's reduce to the ordinary DIS parton distributions. In particular, 

H'^{x,0,0)^q{x), (44) 

where q{x) is the (unpolarized) quark distribution of flavor q. In this limit, 
the Melosh rotation matrices combine to the identity matrix and the following 
simple expression is obtained for the parton distribution 

q{x)=p^Sr^rjfldh6(j2k}j -5(^^-1^) |vI/W({fc,;A„rJ)p. (45) 

This expression agrees with that given, e.g., in Refs. ^ and automatically 
fulfills the support condition, vanishing outside the support region < a; < 1. 
It also satisfies the particle number sum rule, 

Jdxq{x)^Ng, (46) 

where Nq is the number of valence quarks of flavor q, as well as the momentum 
sum rule 

J dxx[u{x) + d{x)] = 1, (47) 
where u{x) and d(x) are the up and down quark distributions. 



5 Results and discussion 

In this Section we present results obtained with two constituent quark models, 
i.e. the relativistic hypercentral quark model of Ref. and the Goldstone- 
boson-exchange (GBE) model of Ref. ]35[| . 

5.1 Adopted models 

The hypercentral model is based on the mass operator M = Mq + V, where Mq 
is the free mass operator. 

Mo = \Jk^+ml (48) 

with ki = 0, and mj being the constituent quark masses. The interaction V 
is taken of the form 

V=-- + Kiy, (49) 
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where y = y + is the radius of the hypersphere in six dimensions and p 
and A are the Jacobi coordinates, 

- ^1 ^ ^2 r ri +r2- 2^3 , , 

The model depends on two parameters, r and k;, and is able to reproduce the 
basic features of the low-lying nucleon spectrum satisfactorily in spite of its 
simplicity. In the hypercentral model the resulting nucleon wave function is a 
product of a space and a spin-isospin part and is SU(6) symmetric. Technical 
details concerning the derivation of the relevant formulae with this model are 
given in the Appendix. 

In the GBE model the same semirelativistic free mass operator Mq of Eq. ( ^8|) 
is combined with an interaction between constituent quarks which is the sum 
of a (linear) color-electric confinement term and a hyperfine (flavor-dependent) 
potential provided by the possible exchange of all mesons of the pseudoscalar 
octet and singlet. Such mesons are considered as Goldstone bosons appearing 
in the model as a consequence of the spontaneous breaking of chiral symmetry. 
The model of Ref. Q depends on five parameters that are fixed by looking at 
the baryon spectrum which is well reproduced up to 2 GeV with the correct 
orderings of the positive- and negative-parity states in the light and strange 
sectors. The resulting nucleon wave functions, without any further parameter, 
yield a remarkably consistent picture of the electroweak form factors . 

5.2 Results 

The spin-averaged {Hi) and the hehcity-flip {Ei) GPD's calculated in the GBE 
model for the u and d flavours are plotted in Figs. ||-^ as a function of x at 
different values of t and They are all positive apart from E'^ that is always 
negative. This is in agreement with the findings within the MIT bag model |Q 
and also within the chiral quark-soliton model |2^, ^ where -t- i?'^ and 
_ £!]d g^j.g leading order in the number of colours Nc, while iJ" — H'^ and 
j^u ^ j^d g^j,g subleading. In all cases the GPD's vanish at x = ^ since in our 
approach they include the contribution of valence quarks only j2^, . They also 
vanish beyond x = \ satisfying the support condition as already mentioned. As 
in the MIT bag model |^ the ^ dependence at fixed t, Figs. ^ and H, turns out 
to be weak, the main effect being a small shift towards larger x with increasing ^. 
In addition an interesting feature emerges from Fig. ^ the large i-independence 
of both and H'' at fixed ^ in the region x > 0.5. Such a result is not consistent 
with the t dependence simply factorizcd in terms of nucleon form factors ]l8| . 
Actually, only the first moments of the GPD's are related to the nucleon elastic 
form factors [Q, i.e. 



j dxH''{x,^,t) ^ F^{t), J dxE'^{x,^,t) ^ F^{t), 



(51) 
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where F^{t) and (^) the contribution of quark q to the Dirac and Pauh 
form factors. This property holds in the present approach, for ^ = only. The ^ 
dependence which emerges from Figs. ^|-^in calculating the first moment of E'^ 
and H'^ is due to the fact that only the valence-quark contribution is taken into 
account. The inclusion of higher order configurations in the Fock space would 
restore the validity of Eq. for ^ > 0. 

It is worthwhile noticing that neglecting the effect of the Melosh rotations in 
the transformation from the canonical- to the front-form wave function, Eq. (j^), 
the helicity-flip GPD's E"^ would be vanishing for wave functions with S waves 
only. The use of the correct transformation can also be appreciated when look- 
ing at the parton distributions for u and d flavours plotted in Fig. ||. Here 
the fully covariant light-front calculation gives contrasting results with those 
obtained with the prescription of Ref. ]30[| including the support and flux factor 
corrections. To preserve normalization the quite substantial shift of the peak to 
lower X values of the light-front calculation is compensated by a higher tail at 
larger x which reflects the presence of high-momentum components in the GBE 
wave functions. 

GPD's calculated with the hypercentral and the GBE model behave quite 
similarly in spite of the fact that the hypercentral model is SU(6) symmetric, 
while the GBE model is not. However, it is known that the SU(6) break- 
ing part in the nucleon GBE wave function is small |3j, In Fig. ^ the 
GPD's obtained within the hypercentral model are shown for ^ = as a func- 
tion of t. In order to better appreciate the difference between the GBE and 
hypercentral CQM's, Fig. ^ shows the t dependence of the ratio between E'' 
and the anomalous magnetic moments k'^ predicted by the respective models 
(A;' = J dx E'^{x,t = 0,^ = 0)). The numerical values for fc' are listed in the 
caption of Fig. ^ and discussed in the text later on. 

In Fig. H the GPD's obtained within the hypercentral model are shown for 
t = —0.5 (GeV)^ as a function of ^. They confirm the results obtained within 
the GBE model. At large x, where it is known that the contribution of va- 
lence quarks is substantial and higher Fock states are less important |2^, the 
distributions are almost independent of ^ so that at fixed t their peak position 
for increasing ^ is shifted to higher values of x by the requirement that they 
have to vanish here at a; = ^. Would higher Fock states (including gluons and 
antiquarks) be included, nonvanishing GPD's should occur for a; < ^ with x = ^ 
as a cross-over point (see, e.g., Ref. p2|). 

Finally we discuss the integral properties of H'' and i?' summarized in 
Eq. (51). In Fig. ^ the proton and neutron Pauli and Dirac form factors are 
shown as function of —t in both the hypercentral and GBE models. The F2 
form factor at i = gives the following values for the nucleon anomalous mag- 
netic moments: k^jgE = 1-20, k^^E = "l-OS, kl^^ = 0.91, k^lyj, = -0.82. 
These results are rather far from the experimental values and consistent with 
analogous light-front calculations when pointlike structure of the quark is as- 
sumed [^3[ The SU(6) breaking effects which produce deviations from 
the value k^/k"" = — 1 come mainly from the Melosh rotations, with a small 
additional contribution from the quark wave functions in the case of the GBE 
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model. The t dependence of Fi and F2 is rather smooth and in both models it 
does not reproduce the experimental results: a clear limitation due to the as- 
sumed pointlike structure of the quarks which can be solved by the introduction 
of specific quark form factors ^ . 

6 Concluding remarks 

A fully covariant approach has been presented linking the overlap representation 
of generalized parton distributions to the non-diagonal one-body density 

matrix in momentum space. As a result of the correct transformation of the 
wave functions from the (canonical) instant-form to the (light-cone) front-form 
description, the support condition is automatically fulfilled and the particle 
number and momentum sum rules are also satisfied. The method has been 
applied to the case of nucleon wave functions involving only the three valence 
quarks in order to study GPD's in terms of constituent quark model wave func- 
tions. This implies that the discussion of GPD's has to be confined to the region 
^ <x < I. In this region one can easily derive the usual parton distributions 
and appreciate the advantages of a fully covariant treatment avoiding problems 
connected with the support condition and the flux factor introduced in other 
approaches. In addition, the helicity flip GPD's naturally arise thanks to the 
relativistic effects of the Melosh rotations. However, the limitation ^ < a; < 1 
prevents the possibility of testing the reduction formula of the first moment of 
the GPD's leading to the nucleon form factors for ^ 7^ 0. 

Results have been presented for the hypercentral model of Ref. and the 
Goldstone-boson-exchange model of Ref. Q . Quite similar results are obtained 
within the two models, the SU(6) breaking contribution present in the GBE 
wave functions being rather small. A strong t and a weak ^ dependence has 
been found in all cases, confirming the results obtained within the MIT bag 
model |2^. Keeping in mind the limitation ^ < x < 1 due to the inclusion 
of the lowest order Fock-space components with three valence quarks only, the 
results obtained show that all the phenomenology for large x and small t can 
be studied within the present approach. 
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A Appendix 

The structure of the nucleon wave function in the hypercentral quark model of 
Ref. 1 27 1 is given by Eq. ( pl| ) where a product is assumed of a symmetric function 



of the momenta with a symmetric function of the spin-isospin variables as in 
Eq. (|4l|). The spin function $f^^(Ai, A2, A3) is symmetric or antisymmetric 
under the interchange of quark 1 and quark 2 for the total spin of the pair 



14 



Si2 = 1 and Si2 ~ 0, respectively. The third quark spin is then coupled to S12 
to obtain the nucleon spin i. The same is done for isospin according to SU(2) 
symmetry. The resulting spin-isospin function, 

'J'Ar(Al,A2, A3,ri,T2,r3) 

= [$°(Ai,A2,A3)$?(Ti,r2,T3) + $l(Ai,A2,A3)$^(ri,T2,T3)] ,(52) 



is then fully symmetric under all permutations. The space part of the nucleon 
wave function is taken with total orbital momentum L — and is written in 
momentum representation as 

V.(fci, fc2, fcs) = V^oo(y)Y[K](f^)' (53) 

where 'ip-y,i^{y) is the hyperradial wave function solution in momentum space of 
an eigenvalue problem for the mass operator M as a sum of the free mass oper- 
ator ( p8| ) and the hypercentral potential given by Eq. (^). The hyperspherical 
harmonics ^[^^j *;^] ^-^e defined on the hypersphere of unit radius. 

Such wave function is transformed to a front-form wave function according 
to Eq. ( p^ ) and inserted into Eq. (|40|). The summation over isospin variables 
gives (5ti20 + ^Ti2i['5r3i/2 + '2S^^_i/2]/3 for the proton and SrriO ^r-i-iji + 

5Ti2i[2(5r3i/2 + '5r3-i/2]/3 for the neutron. 

Summation over spin variables requires much more effort because of the pres- 
ence of the Melosh rotations but is straightforward along the lines of Ref . |Q . 
For the proton we can finally rewrite Eq. (Eg) as 



i=l \i=l ) 

'Kr, {^5?a(« , ^) '5r3l/2 + \X]^S, mr,Xl2 + 25,3 _ 1/2] }(^4) 



where 



Vi({a;J,{Kj_,,}) = 



1 L0\L02'^Z 



-0(fci, fe, fca), (55) 



Mo XxX2X-i 

^va(^ = n^"'(^^)^"'(^')(^'l^3 + *S3-a|A)(Ai^2 + Si -52), (56) 
^"A(fc,fc) = Y^N~\k,)N-\k,)\Y^y{AxA2-Bx-B2)i,y 
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fe 



k k 

B3 J + a, B3^j,-a- B^Sjj, ) | A) , (57) 

with 

N{k) = [(m + xMof + (58) 

= (m + a;-Mo)(m + x^Mo) + /c^/c^ + k'^k^, (59) 

S,,:, = -(m + x[M^)ky + (m + x,Mo)k'y, (60) 

= (m + x'iA'Qk^ - (to + XiMo)k'^, (61) 

B,,,^k'^ky-k'yk^. (62) 

In the above equations the primed (unprimed) kinematical variables refer to the 
final (initial) nucleon state (see Eqs. ( |l0|) and (|lj)) and: 



^o = EV^?+"^'' Af^ = 5]V^'f + mf. (63) 

i i 

In the GBE model the nucleon wave functions are expanded on a basis 
where the spin-isospin part is combined with a space part in the form of cor- 
related gaussian functions of the Jacobi coordinates referring to a particular 
partition. The total wave function is a symmetrized linear combination of such 
basis functions over the three possible partitions thus ultimately violating SU(6) 
symmetry. Specifically, we have 

({fcj, {A.}, {r.}) = c„*„„({fcj, {A,}, {n}), (64) 

n 

where each is summed over all possible partitions, i.e. 

= ^'„(l, 23) + ^-^(2, 31) + *„(3, 12). (65) 
In turn, for each partition (k,pq) 

'^a{k,pq) = ■lpLM{Xk,yk)<^^\{>^pAqAk) j^r{Tp,Tq,Tk), (66) 

where 

i^LM{xk,yk) = xl''+^ yl"+' e-^^'--'y'-+''^~''-y'' Y^Y'ixk, Vk), (67) 
with 2v + X + 2n + I — 2N + L, N being the principal quantum number, and 

Y^i^'ixk, Vk) = Y.^\l^xm\LM)Y^^{xk)Yl^{yk)■ (68) 
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In the above equations Xk,yk are the Jacobi coordinates of partition k. The 
parameters /3, S and 7 are fixed by solving the eigenvalue problem according 
to, e.g., the stochastic variational method as in Ref. |3^. For the nucleon, 

N ^L = 0. 

The calculation with the GBE wave function requires repeating the same 
steps as with the hypercentral wave functions for each partition of the partial 
contribution to the total initial (final) nucleon wave function (|6^). 
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Figure 2: The spin-averaged left panels) and the helicity-flip {E'', right panels) 
generalized parton distributions calculated in the GBE model for flavours u (upper 
panels) and d (lower panels), at ^ = and different values of t: t = (solid curves), 
t = -0.2 (GeV)^ (dashed curves), t = -0.5 (GeV)^ (dotted curves). 
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Figure 3: The same as in Fig. 2, but for fixed t = —0.2 (GeV)^ and different values 
of ^: ^ = (solid curves), ^ = 0.1 (dashed curves), ^ = 0.2 (dotted curves). 
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Figure 4: The same as in Fig. 2, but for fixed t — —0.5 (GeV)^ and different values 
of ^: ^ = (solid curves), ^ — 0.1 (dashed curves), ^ — 0.2 (dotted curves). 




Figure 5: Parton distributions calculated in the GBE model for fiavours u (left panel) 
and d (right panel) within the light-cone approach (solid curves) and the approach of 
Ref. ||3C|] with the inclusion of the support and flux factor corrections (dashed curves) . 
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Figure 7: The ratio of the helicity-flip generalized parton distributions E'' to the 
anomalous magnetic moment fc' calculated in the GBE model (left panels) and in the 
hypercentral model (right panels) for flavours u (upper panels) and d (lower panels), at 
^ = and different values of t = (solid curves), t = —0.2 (GeV)^ (dashed curves), 
t = —0.5 (GeV)^ (solid curves). The values of the quaxk anomalous magnetic moments 
predicted from the two models are: /cqb^ = 1.33, /cgg^ = —0.93, A;^j,p = 1.00, 
kly^ = -0.74. 
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Figure 8: The same as in Fig. m but for the hypercentral model. 
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Figure 9: The Pauli (left panels) and Dirac (right panels) form factors calculated in 
the hypercentral model (dashed curves) and in the GBE model (solid curves) for the 
proton (upper panels) and the neutron (lower panels). 
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